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Abstract: Recent progress in the description of the properties of hadronic atoms on the basis of 
non-relativistic effective Lagrangian approach and Chiral Perturbation Theory (ChPT) is reported. 
For the case of the tt+tt" atom decay, the problem is completely solved, both conceptually and nu- 
merically. For the 7r~p atom, a general expression for the ground-state energy is obtained in the first 
non-leading order in isospin breaking, and the numerical analysis is carried out at 0{p^) in ChPT. 
We briefly consider a possible solution of the "potential model puzzle" in the hadronic atom problem, 
providing a constructive algorithm for the derivation of the isospin-breaking part of the short-range 
hadronic potential from field theory. 



Below, we shall report on the recent progress 
achieved in the description of the hadronic bound 
states - so-called hadronic atoms - within the 
framework of field theory. The reason why such 
a report is delivered at the workshop on heavy 
quark physics, lies in a close similarity of the 
methods that are employed in these, otherwise 
very distinct fields. In fact, it turns out that the 
non-relativistic effective Lagrangian technique that 
was proposed originally by Caswell and Lepa- 
ge 10] to study QED bound states in general, and 
that is widely used now in the physics of heavy 
quarks, provides the most elegant and econom- 
ical approach to the solution of hadronic atom 
problem. In a latter case, one is faced with a 
theoretical challenge of producing a merger of 
the effective non-relativistic theory that implies 
the expansion in the inverse masses of hadrons 
(including Goldstone bosons) and, therefore, in 
the inverse powers of the quark masses, and the 
ChPT that is based on the expansion in quark 
masses. Such a merger can indeed be worked 
out, and the systematic chiral expansion of the 
observables of the hadronic atoms can be ensured 
in this theory. 



* Review of recent work done in collaboration with 
A. Gall, J. Gasser, M.A. Ivanov, E. Lipartia, and 
V.E. Lyubovitskij . 



Recent years have seen a growing interest in 
the study of hadronic atoms. At CERN, the 
DIRAC collaboration |E| aims to measure the 
TT+TT" atom lifetime to 10% accuracy. This would 
allow one to determine the difference oq — 0,2 of 
TTTT scattering lengths with 5% precision. This 
measurement provides a crucial test for the large 
vs small condensate scenario in QCD: should it 
turn out that the quantity ao — 02 is different 
from the value predicted in standard ChPT 0], 
one has to conclude |j] that spontaneous chiral 
symmetry breaking in QCD proceeds differently 
from the widely accepted picture. In the experi- 
ment performed at PSI P, l6|, one has measured 
the strong energy-level shift and the total decay 
width of the Is state of pionic hydrogen, as well 
as the Is shift of pionic deuterium. These mea- 
surements yield isospin symmetric tiN scattering 
lengths to an accuracy which is unique for hadron 
physics. A new experiment on pionic hydrogen 
at PSI has recently been approved. It will allow 
one to measure the decay A^-p -^ ir^n to much 
higher accuracy and thus enable one, in princi- 
ple, to determine the TriV scattering lengths from 
data on pionic hydrogen alone. This might vastly 
reduce the model-dependent uncertainties that 
come from the analysis of the three-body prob- 
lem in Ajr-rf. Finally, the DEAR collaboration M 
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at the DA(f>NE facility plans to measure the en- 
ergy level shift and lifetime of the Is state in 
K~p and K~ d atoms - with considerably higher 
precision than in the previous experiment carried 
out at KEK 1^ for K~p atoms. It is expected (Q 
that this will result in a precise determination of 
the / = 0, 1 S'-wave scattering lengths. It will be 
a challenge for theorists to extract from this new 
information on the KN amplitude at threshold 
a more precise value of e.g. the isoscalar kaon- 
sigma term and of the strangeness content of the 
nucleon. 

In order to carry out the precision experi- 
mental tests of QCD mentioned above, on the 
theoretical side one faces the problem of finding 
the relation between the measured characteris- 
tics of hadronic atoms - energy levels and decay 
probabilities - and the strong hadronic scattering 
lengths in the isospin limit. In general, we have 
the following relations between these quantities 

A£;,tr-*^Reacc(l + 5.), (1) 

Fco ~ (phase space) x ^1\aco{l + 5r)\'^ ■ 

Here Ai<^str denotes the strong energy-level shift 
in the Is state from its Coulomb value (total shift 
minus the pure QED contribution) , and Vco is the 
partial decay width into the neutral channel (e.g. 
A^-p — > Tr'^n); Ucc and Oco stand for the par- 
ticular isospin combinations of the strong scat- 
tering lengths (charged and neutral channels are 
marked by subscripts "c" and "0", respectively), 
and ^0 denotes the value of the Coulomb wave 
function at the origin. The quantities S^ and 6r 
stand for the isospin-breaking corrections and, in 
general, depend on the details of the strong dy- 
namics. In order to extract the scattering lengths 
from the experiments, these quantities must be 
known to an accuracy that matches the accuracy 
of the measurements. 

Historically, the potential model was the first 
one been applied to the calculation of the had- 
ronic atom characteristics. In brief, this model 
assumes that the Coulomb effects and the mass 
differences between the charged and neutral par- 
ticles in the same multiplet are responsible for all 
isospin-breaking effects in the low-energy hadron 
physics: the short-range potential that describes 
strong interactions, is presumed to be isospin- 



symmetric. In this manner, one may calculate 
the isospin breaking corrections 5^ and Sp |5|, 
g|. The results of these calculations are how- 
ever in the striking disagreement with the field- 
theoretical evaluations based on ChPT (see be- 
low) . It turns out that the assumption about the 
isospin symmetry of the short-range hadronic po- 
tential is too restrictive - one should allow for a 
small isospin-breaking piece in the short-range 
interactions as well, in order to achieve an agree- 
ment with the field-theoretical calculations. In 
this way - providing a constructive algorithm for 
the derivation of the isospin-breaking piece from 
field theory - one may achieve the solution of ap- 
parent potential model puzzle. 

Recently, using a non-relativistic effective La- 
grangian framework, a general expression for the 
decay width Tj^^^^t^Oj^o of the Is state of the 
TT+TT" atom was obtained at next-to- leading or- 
der in isospin-breaking mn. Numerical analysis 
of this quantity was carried out at order 0{e^p^) 
in ChPT |11 . These investigations have con- 



firmed and generalized the results of earlier stud- 
ies ||lj, ^. The expression for the decay width 
has the form 



-4^7r = ao - a2 + 0[a, {rud ~ rriu)'^) ., 



K^ 



AAf2 



2a 



^ (ao + 202)^ (In a - 1) X 



9Af^ V " ■ -' 3 



(2) 



X (2ao -I- 02) + o{a, {md - to„)^) . 

Here p* = {M^+ - M^,, - ^M^^a^^^, and a/, 
(/ — 0, 2) denote the strong tttt scattering lengths 
in the channel with total isospin /. The quan- 
tity Atttt is calculated as follows pffl. One cal- 
culates the relativistic amplitude for the process 
TT+TT^ -^ TT^TT^ at 0(a, (md — ruu)'^) in the nor- 
malization chosen so that at 0(1) the amplitude 
at threshold coincides with the difference oq — 02 
of (dimensionless) 5'-wave tttt scattering lengths. 
Due to the presence of virtual photons, the am- 
plitude is multiplied by an overall Coulomb phase 
9c that is removed. The real part of the remain- 
der contains terms that diverge like |p|~^ and 
In2|p|/Af^+ at |p| ^0 (p denotes the relative 3- 
momentum of charged pion pairs). The quantity 
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Atttt is obtained by subtracting these divergent 
pieces, and by then evaluating the remainder at 
p = 0. 



Re(e-^^=i,,) 



IpI 



boM 



2|P| 



M^ 



Stt 



■ J\-n 



(3) 



As it is seen expUcitly from Eq. (^, one can 
directly extract the value of A-n-n from the mea- 
surement of the decay width, because the cor- 
rection Kt^t^ is very small and the error intro- 
duced by it is negligible. We emphasize that in 
derivation of Eq. (0), chiral expansions have not 
been used. On the other hand, if one further 
aims to extract strong scattering lengths from 
data, one may invoke ChPT and to relate the 
quantities A-^t^ and oq — 02 order by order in the 
chiral expansion. This requires the evaluation 
of isospin-breaking corrections to the scattering 
amplitude. At order O(e^p^) in chiral expansion 
we obtain [Q (for the values of scattering lengths 
ao = 0.206, 02 = -0.0443) 



A^^=ao-a2 + e, e = (0.58 ± 0.16) • 10"^ 
^ 5t = 0.058 . (4) 



K = 1.07- 10" 



At this stage, one may recall that the calcula- 
tions in the potential model with the isospin- 
symmetric strong potential yields S-p being of neg- 
ative sign and the same order of magnitude 0] . 
In the case of the 7r~p atom, the treatment 
proceeds along the lines very similar to those for 
TT+TT" case [|lj]. Our investigations are aimed at 
the derivation of the general expression for the 
■K^p atom energy-level shift in the Is state. The 
total shift is given by a sum of the electromag- 
netic (pure QED) and strong pieces. Our cal- 
culations for the electromagnetic shift within a 
high accuracy yield the same result as given in 
Ref. [Q. The final result for the strong shift in 
the first non-leading order in isospin breaking is 
given in a form similar to Eq. (Eh 



AE^str = -2a>^ A^N (1 + K^n) 



(5) 



Kt^IsI = 2a/ic(l — In a)A.nN + o(a, rud — ruu) , 

where ^c denotes the reduced mass of the 7r~p 
pair, and the quantity A^n is defined analogously 



to Atttt- To calculate this quantity, one has to 
evaluate the 7r~p — > ■K~p relativistic scattering 
amplitude at 0{a,md — rriu), drop all diagrams 
that are made disconnected by cutting one pho- 
ton line, and discard the spin-flip piece. The re- 
mainder is denoted by ij^N- The regular part 
of Ittn at threshold defines the quantity AjtN in 
analogy to Eq. (||) 



Re ie-'^'^ Un) 



IpI 



B2 In ■ 



Mc 



27r 



•^TT. 



N, 



(6) 



and the normalization of the amplitude is chosen 
so that Attn = bo — bi + 0{a,md — my), where 
60 and 61 denote the isospin even and odd strong 
ttN scattering lengths. 

In order to extract the value of &o ~ ^1 from 
the Ti^p measurement, one may again resort to 
ChPT, to calculate the isospin-breaking correc- 
tions to the 'kN scattering amplitude at thresh- 
old. At chiral order 0{p^) where only the tree 
diagrams contribute, the result looks as follows 



A. 



(2) 



■N 



&0 - &1 + e 



(2) 
•Af 



.(2) 
-■kN 



mp(8ciA^ - 4eVi - 6^2) 
87r(mp-HM^+)F2 



(7) 



where A^ = M^^ — M^o denotes the pion mass 
difference, and ci (fi) are the strong (electro- 
magnetic) low-energy constants (LEC's) from the 
0(p2) Lagrangian of ChPT ||l^. In order to per- 
form the numerical analysis, one has to specify 
the values of these LEC's. The "strong" constant 
ci can be determined from the fit of the elastic 
ttN scattering amplitude at threshold to KA86 
data fill: a ^ -0.925 GeV"^ The value of the 
constant /2 can be extracted from the proton- 



neutron electromagnetic mass difference [17 



6^2 == (-0.76 ± 0.3) MeV. The determination 
of the constant /i from data is however, prob- 
lematic. For this reason, in our analysis we have 
used order-of-magnitude estimate for this con- 
stant: — I/2I < /i < 1/2 1- With these values of 
the low-energy constants, we obtain the isospin- 
breaking correction to be S^ = (— 4.7±2.0) • 10~^, 
again in striking disagreement from the potential 
model prediction ^6, = (-2.1 ± 0.5) • lO^^. It 
remains to be seen, how the 0{p^) results are 
altered by the loop corrections in ChPT mM. 
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Given a systematic discrepancy of the po- 
tential model predictions with the results of cal- 
culations based on ChPT, it is natural to seek 
a derivation of the potentials that are used in 
the potential model, on the basis of ChPT. In 
a slightly more restricted context, one may ask, 
how the isospin-breaking part of the short-range 
"strong" potential is obtained from ChPT, when 
the isospin-symmetric part is already known to 
fit well ChPT predictions (we recall that the iso- 
spin-breaking part is assumed to vanish identi- 
cally in existing potential models |^, |^). It is 
widely presumed that the potential constructed 
from the field theory will be necessarily singular 
in the position space and will require some kind 
of regularization. Based on a simple solvable 
model, we shall however demonstrate that this 
is not the case: almost any well-behaved short- 
range potential, including those that were used 
in Refs. g |^, can be generahzed to include prop- 
erly the full content of isospin-breaking effects in 
ChPT. 

The key to the solution of the problem given 
above, lies in the universality conjecture. This 
conjecture - completely in spirit of the low-energy 
effective Lagrangian approach to bound systems 
- states that the bound-state energies in the field 
theory, and in the potential model are the same 
at the first order in isospin breaking, provided 
the quantities A-„n calculated in these two the- 
ories, coincide. We shall ensure the universality 
for the case of a simple model where the interac- 
tion Hamiltonian is given by a sum of Coulomb 
and short-range rank-1 separable interactions 

C/(p, k) = C + V = - j^^ -f Az;(p)z;(k) , 



|p-k|^ 



u(q) 



/3^ 



/32-Hq2 



(8) 



where (3 denotes the cutoff mass. The general- 
ization to the case of generic potentials, multi- 
channel case, inclusion of relativistic effects, is in 
progress and will be reported elsewhere Q . 

With a given interaction potential, one may 
evaluate the energy- level shift of the ground state 
of the bound system. The equation for the po- 
sition of the bound-state pole in the (complex) 



energy plane is given by [[lO[ 

z-.Bo-(*o|r(z)|*o)-0, 



(9) 



t{z) - V + V(z - Ho)-i(l - |vI'o)(vI'o|)r(z) . 

It is straightforward to solve this equation in the 
perturbation theory. The solution at the first 
non-leading order in a reads 

4a/i 



Ai? = *^ro(l-^ 



-t- To [in a 

IT 



-4l 



To = A 1 



mA/3 
47r 



(10) 



and /i stands for the reduced mass. 

Next, we calculate the regular part of the 
elastic scattering amplitude at threshold, with 
the normalization according to Eq. (|6|) 



A 



^Inl^ (11) 



From Eqs. (O) and (11) one immediately ob- 
tains that the energy level shift in the potential 
model is given again by Eq. (ra) - that is, the uni- 
versality holds in that particular case, considered 
here. 

Based on the universality conjecture, we can 
provide a constructive algorithm for the deriva- 
tion of the isospin-breaking part of the short- 
range potential V from ChPT. The amplitude 
at threshold in the latter is generally given by 



A = Aq + Ai 



where -4o(i) denote the 



isospin-conserving(breaking) parts of the ampli- 
tude, and ellipses stand for higher-order terms in 
isospin breaking. In order to ensure the inclusion 
of the full content of isospin-symmetry breaking 
in ChPT into the potential model, it thus suf- 
fices to match the amplitude A in both theories. 
The problem evidently has too much freedom, 
and we can choose one parameter in the poten- 
tial - say, the coupling constant A - to obey the 
matching condition for the amplitude. Writing 
A = Ao + Ai + ■ • •, the matching condition yields 

27r _^o 



Ao 



Ai = - 



1 + 
2tt Ai 



■^0 



(12) 



2^iaAl In ^ 



Uo)2 
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The matching condition (|l2|) solves our prob- 



[13] V.E. Lyubovitskij and A.G. Rusetsky, Phys 



lem completely - the bound-state energies calcu- 



Lett. B 389 (1996) 181; V.E. Lyubovitskij, 
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coincide, by definition, with those calculated on 
the basis of ChPT. Wo hope that — after the suit — 
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